Metal-Insulator phase diagram and orbital selectivity in 3-orbital models with 

rotationally invariant Hund coupling 



OO 
O 

o 

O 

Q 

00 



I 

C 

o 
o 



> 
o 

in 

(N 
00 

o 

> 

X 



Philipp Werner, 1 Emanuel Gull, 1 and Andrew J. Millis 2 

1 Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland 
2 Columbia University, 538 West, 120th Street, New York, NY 10027, USA 
(Dated: December 8, 2008) 

A three band model containing the essential physics of transition metal oxides with partially filled 
t2 g shells is solved in the single-site dynamical mean field approximation, using the full rotationally 
invariant Slater-Kanamori interactions. We compute the metal-Mott insulator phase diagram in the 
space of chemical potential and interaction strength, determine the response of the different phases 
to perturbations which break the orbital symmetry, and establish the regimes in which an orbital 
selective Mott phase occurs. The results are compared to data on titanates, ruthenates, vanadates 
and Ceo- 
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I. INTRODUCTION 

The physics of strong (electronic) correlations plays a 
central role in modern-day condensed matter physics. 1,2 
The essence of this problem is the competition between 
the repulsive interactions felt by electrons in transition 
metal d-orbitals or lanthanide/actinide / orbitals and 
the itineracy arising from hybridization with other or- 
bitals in the material. For an atom in free space, the 
d and / shells have respectively a 5 and 7-fold or- 
bital degeneracy and when the orbitals are partially 
filled, Coulomb interaction effects lead to a compli- 
cated multiplct structure. In a solid state environment 
the orbital degeneracy may be fully or partially lifted. 
In some cases, for example the cuprate high temper- 
ature superconductors^ the degeneracy is fully lifted 
and the low energy physics may be described by the 
one-band Hubbard model in which the multiplet struc- 
ture is triviali^ However, for many materials of inter- 
est, including for example the (La/Ca)Ti03 series, the 
SrVOa; materials, the (Sr/Ca)Ru03 compounds and their 
Ruddlesden-Popper variants, the new Fe based supercon- 
ductors and the A„C6o series of materials, as well as es- 
sentially all interesting lanthanide/actinide compounds, 
the orbital degeneracy is not fully lifted and nontrivial 
multiplct effects are expected to be important. Of par- 
ticular interest is the effect of orbital degeneracy on the 
Mott metal-insulator transition. It is generally believed^ 
that the critical interaction strength required to drive 
a metal-insulator transition depends on the orbital de- 
generacy, being larger for systems with several degener- 
ate orbitals than it is for one-orbital models. This gives 
rise to the physics of orbital selectivity, whereby a bro- 
ken orbital symmetry, either spontaneous or induced by 
crystal symmetry, may drive some orbitals into insulating 
states. This phenomenon has for example been argued to 
be of crucial importance in understanding the insulating 
phase of LiTiOyi and of the metal insulator transition in 
Ca 2 Ru0 4 ^ 

The dynamical mean field theory (DMFT) provides a 



non-perturbative method to study the interplay of cor- 
relation effects and electron banding and has in particu- 
lar produced insights into the correlation-driven (Mott) 
metal-insulator transition in the one-orbital models 
While the issue of the Mott transition in multiorbital 
systems has been addressed by various techniques ) 10 ' 11 : 12 
a comprehensive picture has not emerged, in part be- 
cause of the theoretical difficulties associated with the 
treatment of the various Hund and pair hopping terms 
required for a realistic treatment of partially filled d- 
orbitals. Dynamical mean field theory maps a lattice 
problem onto a quantum impurity model (a finite size 
system coupled to a noninteracting bath of electrons) 
plus a self consistency condition. For systems in which 
the orbital degeneracy is fully lifted the quantum im- 
purity model is a variant of the one-orbital "Anderson 
Impurity Model" , for which powerful numerical tech- 
niques have been known for many years i 13 ' 14 ' 15 : 16 : 17 How- 
ever, these techniques encounter difficulties when ap- 
plied to materials with partially filled, degenerate d- 
orbitals, where the on-site interaction includes both spin 
exchange and "pair hopping" terms. The Hirsch-Fye 
method, which has been the standard approach for multi- 
orbital models with density-density interactions, relies on 
a Hubbard-Stratonovich transformation of the interac- 
tion term. In the orbitally degenerate case the multiplic- 
ity of interactions requires many auxiliary fields, which 
become difficult to sample. Rotational invariance be- 
comes very difficult to preserve and a severe sign prob- 
lem is reported^ The proliferation of states also creates 
difficulties for exact diagonalization methods, although 
recent progress has been made along these line . 11 : 12 

In this paper we exploit a recently develope d 19 : 2 im- 
purity solver which is free from the defects of the other 
methods. In this method the on-site Hamiltonian is 
solved exactly, and the coupling to the bath is treated by 
a perturbation expansion which is sampled stochastically 
via an importance-sampling procedure. The method al- 
lows a detailed and accurate treatment of thermodynamic 
quantities and (via analytical continuation) of dynamics, 
down to temperatures of the order of 0.1% of the basic 
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energy scales of the problem. An additional benefit of 
the method is that it provides information about which 
configurations of the correlated site make the dominant 
contributions to the partition function. While the com- 
putational effort of our method scales exponentially with 
the number of orbitals, it can easily handle three orbitals 
on desktop machines, and five orbitals on larger clusters. 

In this paper we use the method to analyze the "three 
orbital" model which is relevant to materials such as 
LaTi03, SrVC>3 and SrRu03, where the physics is dom- 
inated by electrons in the transition metal ti g orbitals. 
The model is also relevant to electron-doped Ceo, where 
the three orbitals correspond to the triplet of LUMO 
states of C 6 o- We determine the metal-insulator phase 
diagram, study the response to perturbations which lift 
the orbital degeneracy and determine the orbital selectiv- 
ity of the doped Mott insulating state. Our work builds 
on our previous investigation of a "two-orbital" model^ 
relevant to systems with e g symmetry. In the two orbital 
case, in the presence of strong Hund coupling, one has ei- 
ther a one-electron state or a filled (spin-polarized) shell. 
The new feature of the three orbital model is the case 
n = 2, where one can have a multiclcctron state with 
high local spin alignment but a partially filled shell. 



II. FORMALISM 

We study a model involving three orbitals (labeled by 
a = 1, 2, 3), with Hamiltonian 

Here i labels sites in a lattice and k a wave vector in the 
first Brillouin zone, ni >a ,a = d\ a a di t a,a is the density of 
electrons of spin a in orbital a on site i, /i is the chemical 
potential, A a is a level shift for orbital a arising from 
a ligand field splitting and e£° is the band dispersion. 
For the following analysis, the relevant property of the 
dispersion is the density of states N(oj) given by (j (dk) 
symbolizes an integral over the appropriate Brillouin zone 
with the correct measure factors) 

N ab (u,) = j (dk)8{u - ef). (2) 

We have assumed that the symmetry is such that the lo- 
cal density of states is orbital-diagonal and independent 
of a; this assumption holds for pseudocubic materials 
such as the La-titanates and the "113" Sr/Ca ruthenates, 
as well as for A„C6o- We expect that the qualitative con- 
sequences of a symmetry breaking in the density of states 
are similar to those obtained by introducing an explicit 
crystal field splitting A a . 

For the interaction term we take the standard Slater- 



Kanamori form (we have suppressed the site index) 

-ffloc = E Un aA n a.,l 
a 

+ ^2 \u' n a,a n b-a + (U' - J)n aj(T n b ^ 

a>fc,tx 

_ E J ( d a,l d b,t d b,l d a,'\ + d b,l; d tl da <1 da 'i + h - c -)- 

(3) 

Here U is the intra-orbital and U' the inter-orbital 
Coulomb interaction, while J is the coefficient of the 
Hund coupling and pair-hopping terms. We adopt the 
conventional choice of parameters, U' = U — 2J, which 
follows from symmetry considerations for d-orbitals in 
free space and is also believed to hold in solids. With this 
choice the Hamiltonian ([3]) is rotationally invariant in or- 
bital space. The chemical potential required to obtain a 
given occupancy at fixed U decreases as J is increased; 
for example the condition for half filling is jj, — |Z7 — 5 J. 
We shall focus on the case U > 3J, in which (loosely 
speaking) the U interaction controls the occupancy and 
once the occupancy is fixed the J interactions then con- 
trol the arrangement of the electrons among orbitals. For 
U<3J the physics is different: the local level first max- 
imizes the spin and then adjusts the local occupancy ac- 
cordingly. We are not aware of materials for which this 
regime is relevant. 

We solve the model using the single-site dynamical 
mean field approximation, 9 which neglects the momen- 
tum dependence of the self-energy and reduces the orig- 
inal lattice problem to the self-consistent solution of a 
quantum impurity model given by the Hamiltonian 

Hqi = - EG" _ &a)n a , a + Hloc + ffhyb + #bath (4) 
a,rr 

with 

#hyb = E V k,a,<r d l,,<T C k,a,,<T + h.C, (5) 
ffbath = E £ a ath ( fc ) C L, CT C fe ,a,.. (6) 

The important quantity for the subsequent analysis is 
the hybridization function A^(w) which depends on or- 
bital a, spin er and frequency and whose imaginary part 
is 

ImA»£M = J{dk)\V h , a , a \ 2 8{u-e h a ^{k)). (7) 

In the computations presented here we take an orbital- 
independent semi-circular density of states with band- 
width At (Bcthe lattice). The t2 g band widths for early- 
stage transition metal oxide compounds are of the order 
of 3eV, so that t - 0.75eV. 

The hybridization function is fixed by a self- 
consistency condition^ involving the impurity model 
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FIG. 1: (color online) Electron density n per orbital per spin 
computed as a function of chemical potential p for different 
values of the interaction parameter J at U/t = 12 and tem- 
perature fit = 50. The orbital symmetry of the Hamiltonian 
is unbroken (A a = 0) and orbital and spin symmetry were 
enforced in the calculation. Plateaux in n(p) correspond to 
Mott insulating states. Open (full) symbols correspond to 
metallic (insulating) solutions. 



Green's function G^ 1 , the self energy E QI of the quan- 
tum impurity model and the momentum integral of the 
Green's function of the lattice problem computed with 

sQ 1 , 



de- 



N a (e) 



iuj n - e - 



(8) 



We note that insulating solutions may easily be distin- 
guished from metallic solution by the behavior of G^ 1 (t): 
for an insulator at low T this quantity drops exponen- 
tially as r is increased from or decreased from /3, while 
in a metallic phase the constant Fermi level density of 
states leads to a slow power-law decay. 

The simulations were performed using a continuous- 
time QMC solver which samples a diagrammatic expan- 
sion of the partition function in powers of the impurity- 
bath hybridization ffhyb' 19,20 We monitored the particle 
densities in each orbital, the Green's functions and self 
energies of the impurity model, and the contributions of 
each eigenstatc of H\ oc to the partition function. For a 
three-orbital model the dimension of the Hilbcrt space of 
H\ oc is 64, so the series is constructed in terms of traces 
of products of 64 x 64 matrices combined with determi- 
nants made up of the the hybridization function Ahyb 
evaluated at different time arguments. The bottleneck of 
the simulation is the trace computation. To speed this 
up, it is important to group the eigenstates of Hi oc ac- 
cordin g to the conserved quantum numbers as explained 
in Ref. |22| . The matrix-representation of the operators d 
and dft then acquires a block structure, with blocks of size 
< 9. This way, the simulation becomes efficient enough 
to run on a desktop machine. Our results were obtained 
using about 3-5 CPU hours per iteration. 



III. METAL-INSULATOR PHASE DIAGRAM; 
ORBITALLY SYMMETRIC CASE 



To map out the metal-insulator phase diagram we have 
computed the dependence of density (typically repre- 
sented as a density per orbital per spin) as a function 
of chemical potential for various interaction values. Fig- 
ure [1] shows representative results. For sufficiently neg- 
ative /z the solution we find has density n = ("band 
insulator"). As fi is increased, the density increases. For 
small U the increase is smooth at all fx, while at larger 
U plateaux occur at which the density is fixed to the in- 
teger values n — 1,2,3 (so the density per orbital per 
spin is fixed to 1/6,2/6,3/6). We identify the regions in 
which n smoothly increases as metallic phases and the 
plateaux as Mott insulating regions; we have confirmed 
these identifications by examination of G^ 1 (t). Metallic 
(insulating) solutions are plotted with open (full) sym- 
bols. At the U = 12t value studied in Fig. [T] we see that 
for J/t — and 1, we have plateaux at each integer n, 
for J/t = 2 only at n — 2, 3 and for J/t = 3 and 4 there 
is only a plateau at n = 3. From similar traces at var- 
ious values of U, J we have constructed metal-insulator 
phase diagrams in the plane of chemical potential fj, and 
correlation strength U. 

The upper panel of Fig. [2] shows the phase diagram for 
J = 0. We see that the critical U required to drive a 
Mott transition depends weakly on density, ranging from 
U = 6.5t at n = 1 to U = lOt at n = 3. Positions and 
widths (in /i) of the Mott lobes are only weakly dependent 
on band filling (at U — 16i the width is about lit for all 
three lobes). The lower panel of Fig. [2] shows that the 
situation changes quite dramatically in the presence of 
a Hund coupling. The size of the 3-electron insulating 
lobe is substantially increased at the expense of the 2- 
and 1-electron lobes. Furthermore, the value of U C 2 for 
the half-filled insulating state is reduced from w lOt to 
w 3t, while the 2- and 1-electron insulating lobes shift to 
higher values of U. 

Insight into the physics of the metal-insulator phase 
boundaries can be obtained by considering the atomic 
limit. If E n denotes the lowest eigenvalue of the n- 
particle sector of H\ oc , then an estimate of the Mott gap 
is 



A 



Mott 



(n) — E n+ i + E n -i — 2E n 



(9) 



The actual Mott gap is reduced by an amount of the 
order of the bandwidth (At) while the critical U required 
to drive a metal-insulator transition may be estimated 
by comparing the strong coupling Amou to the electronic 
kinetic energy K = — £ a (T j {dk)e k (c{ a a c k ^ a ) . 

For J = the interaction term H\ oc may be rewritten 
H loc = UN tot (N tot - l)/2, so Eq. © gives A Mo tt = V 
for all n. The upper panel of Fig. [2] shows that at large 
U the width in \x is almost the same for the three Mott 
lobes, consistent with this simple argument. Similarly 
the n-dependence of the critical U is consistent with the 
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FIG. 2: Metal-insulator phase diagram presented in the space 
of chemical potential fi and interaction strength U (measured 
in units of the quarter-bandwidth t) for A a = 0, fit = 50 at 
Hund's coupling J — (upper panel) and J = U/6 (lower 
panel). Orbital and spin symmetry were enforced in the cal- 
culation. Error bars are of the order of the symbol size. The 
numerals in the lobes indicate the electron concentration per 
site in the insulating phases. In the lower panel the solid di- 
amonds indicate the boundary of a spin freezing transition 
discussed in Ref. |28| . while the line with squares plots the 
locus of fi and U corresponding to the density n = 2. 
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FIG. 3: Phase diagram in the plane of Slater-Kanamori pa- 
rameters U and J calculated for the orbitally symmetric 
model (A a = 0) at fit = 50 and half filling (n = 3). The 
hashed region (U < 3J) corresponds to an effectively attrac- 
tive Coulomb interaction; this situation does not normally 
occur in transition metal oxides and is not considered here. 



computed for fully spin polarized electrons, and should 
be reduced by a factor of y2 to account for the double- 
exchange physics of spin polarized electrons hopping in 
a paramagnetic environment. These considerations give 
K(n = 1) « K(n = 2) w 0M; the reduced K and 
reduced U account for the shift of the critical U. For 
n = 3, the situation is different: as J becomes large, the 
constraint of total on-site spin polarization means that 
no low-energy states are available for conduction: there 
is only virtual hopping and as in the half-filled double ex- 
change model one would have insulating behavior driven 
by J only. This means that the kinetic energy is very 
rapidly suppressed by a non-vanishing J, explaining the 
rapid shift in the phase boundary seen in Fig. [5] and in 
more detail in Fig. O 



n dependence of the noninteracting kinetic energy K(n = 
1) w lAt, K(n = 2) w 2.3f and K(n = 3) w 2.5*. 

A non-zero J term changes the energetics. The lowest 
energy state is of maximal spin and we find Ajviott(l) — 
Amou(2) = U — 3 J. However, for n = 3, adding a fourth 
electron requires flipping a spin, so that Amo«(3) = 
U + 4 J. For J = U/Q as in the lower panel of Fig. [5] this 
becomes A Mo tt(« = 1,2) « U/2 and A Mo tt(3) ~ 5Z7/3. 
These considerations explain the comparable widths of 
the Mott lobes for n = 1,2 and the much larger width 
of the n = 3 Mott lobe. The variation of the critical 
U is more subtle. For n = 1,2 the transition occurs 
at a sufficiently strong correlation that we may assume 
that each site is always in its maximal spin state, al- 
though our calculation is in the spin-disordered phase, 
so the direction of the moments is random from site to 
site. The noninteracting kinetic energy should then be 



IV. LIFTING OF THE ORBITAL 
DEGENERACY: METAL-INSULATOR PHASE 
DIAGRAM AND MOTT INSULATING STATES 

In this section we consider the consequences of an ex- 
plicit breaking of the orbital symmetry of the model. We 
focus mainly on J > and dopings between n = 1 and 
ri = 3. The cases of n = 1 and n = 3 are straightforward. 
At J > 0, the n = 3 state is a filled shell, stable against 
orbital splitting for small differences among the A Q , while 
for larger crystal field splitting a high-spin/low-spin tran- 
sition will occur, with physics analogous to that discussed 
in the two orbital context in Ref. j21|. For n = 1 the 
qualitative behavior is clear: the model becomes cither 
an effective one orbital model or an effective two orbital 
model; the physics of these two cases has been previ- 
ously discussed^ 1 As an example we show in Fig. H] 
the evolution of the orbital occupancy under a "trigo- 
nal" crystal field splitting A a = A, A 2 = 0, A 3 = -A 
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FIG. 4: Orbital filling as a function of crystal field splitting 
A in the symmetric case: Ai = —A, A2 = 0, A3 = A. 
The parameters are U/t = 7 < U c2 , J/U = 1/6, fit = 50, 
and the density at A = corresponds to 1 electron. As the 
crystal field splitting is increased, band 1 (which is raised) 
empties out and undergoes a metal-band insulator transition 
near A/t ~ 0.4. At the higher value A « 0.8t the second band 
empties out, leaving what is effectively a one orbital model 
for which U > U C 2 so the state is insulating. 
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FIG. 5: Effect of a "1 up, 2 down" crystal field splitting on 
the 2-electron insulating phase. Heavy black line: boundary 
of the two electron Mott insulating state in the space of in- 
teraction U and chemical potential fi computed for Ai = t, 
J = [7/6 and fit — 50. The crystal field splits the threefold 
degenerate level into a doublet and a singlet, with the singlet 
lying higher. Dashed lines: metal-insulator phase boundary 
for the same interaction parameters and A a = for compar- 
ison. 



which separates all three orbitals. We see that as the 
splitting is increased one band becomes depopulated and 
then, at a higher A, the second band empties out, leaving 
an orbitally polarized Mott insulator. This behavior is 
consistent with the proposal of Pavarini et aL,— who ar- 
gued that the insulating behavior of the n = 1 material 
LaTiC>3 is due to a relatively strong ligand field which 
splits the degeneracy of the three orbitals. 



We focus now on the case n = 2, which is relevant for 
example to SrRu03 and the Ruddlesden-Popper materi- 
als Sr^+iRunOsn+i. In studying these cases our main 
focus is on the simplest symmetry breaking, a cubic- 
tetragonal distortion which splits the 3-fold degeneracy of 
the tig state into a singlet and a doublet. We parametrize 
this splitting by moving one orbital (which we take to 
be "orbital 1" by an energy A while keeping the other 
two fixed, so Ai = A, A 2 = A 3 = 0. There are two 
cases: either the doublet lies lower than the singlet ("1 
up, 2 down", A > 0) or the reverse ("1 down, 2 up", 
A < 0). While we implement here the symmetry break- 
ing by shifting the orbital energies, other ways of break- 
ing the symmetry (e.g. choosing different bandwidths) 
will have similar effects. 

We begin by considering the large-J7, Mott insulating 
regime. Figure [5] compares the metal insulator phase 
boundary computed for the orbitally symmetric model 
to the location of the n = 1, 2 Mott lobes computed for a 
"1 up, 2 down" crystal field Ai of magnitude t. Magnetic 
and orbital ordering are suppressed. Lifting the orbital 
degeneracy is seen to have a very substantial effect on 
the n = 2 Mott phase and a noticeable but less dramatic 
effect on the n = 1 Mott lobe. The critical interaction 
strength needed to drive the two electron phase insulat- 
ing is seen to be reduced to less than half of the value 
found in the orbitally symmetric model. The width of 
the 2-electron insulating plateau is enhanced, but to a 
lesser extent: the increase in the width is approximately 
Ai. Both positive (1 up, 2 down) and negative (1 down, 
2 up) crystal field splittings stabilize the insulator, but 
the effect of a positive Ai (which shifts band 1 up) is 
much larger. For Ai = — t (not shown) the end point of 
the 2-electron lobe is U C 2 ~ 9.7t. The difference occurs 
because if one level is shifted up, the n = 2 electron state 
effectively becomes a filled shell which (as can be seen for 
the 3 electron state in Fig. ^ is particularly stable. 

Figure [5] presents the response of the 2-electron insu- 
lating state to crystal field splitting for two values of J. 
In the high-spin filled shell case of two electrons in two or- 
bitals studied in Ref. [2l], the insulating state (for J > 0) 
did not respond at all to a weak crystal field splitting. 
Here, because at density n = 2 the A a = state is not 
a filled shell, the two electron insulating state responds 
even to an infinitesimal crystal field splitting: the orbital 
susceptibility is non- vanishing. 

The value of the Hund coupling J has important effects 
on the response to a crystal field. At J = (rightmost 
traces in Fig. [5] there is no energetic barrier to placing 
two electrons in the same orbital. If Ai is decreased 
(/i — Ai increased) the occupancy of band 1 increases to 
1 per spin while the occupancy of the other two bands 
decreases smoothly to zero. If Ai is increased, the non- 
degenerate state empties out while the occupancy of the 
two degenerate states remains equal, and approaches 1/2 
per spin per orbital. If J > (left hand traces) the sit- 
uation changes: at A a = the lowest energy state is 
spin triplet, so that as Ai is decreased only 1/2 electron 
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FIG. 6: (color online) Orbital filling as function of crystal field 
splitting Ai computed for the two electron insulating state 
with U/t — 12, /3t = 50 and indicated values of J/t. In order 
to display all of the curves on the same figure the crystal field 
coordinate is chosen to be p — Ai where the value of jj, corre- 
sponds to n/t=16, 14 for J/t=0, 1, respectively. The ground 
state is insulating for all points shown. The red lines with 
circles correspond to the occupancy of orbital 1 and the blue 
lines with stars to the occupancy of orbitals 2, 3. The conven- 
tions are such that increasing Ai to positive values (moving 
to the left on the plot) shifts the non-degenerate orbital up 
(decreasing its occupancy). The offset between the curves for 
different J values arises because of the J-dependence of the 
location of the Mott lobes. Dashed horizontal lines are shown 
at filling 1/4 and 1/2. 



per spin can populate orbital 1, which leads to an aver- 
age distribution of (1/2, 1/4, 1/4). For even larger Ai a 
high-spin/low spin transition will occur, but we do not 
consider this here. Depending on the degree of band- 
nesting, the states considered here may become unstable 
to orbital ordering. Magnetic and orbital ordering will 
be discussed in a future paper ^ 



V. CRYSTAL FIELDS AND THE DOPED MOTT 
INSULATOR 

We now consider the behavior occurring as the n = 2 
Mott insulator is doped in the presence of a non- vanishing 
crystal field splitting. The results presented in this sec- 
tion pertain to an orbitally symmetric solution, and may 
be changed if orbital order occurs. Our preliminary re- 
sults are that except very close to the Mott insulating 
phase boundaries, the doped states are stable against 
staggered orbital ordering^ 

Representative data are shown in Fig. [7] which plots 
the dependence of orbital occupancy on chemical poten- 
tial for a relatively small (|Ai| = 0.25i, upper panel) and 
relatively large (|Ai| = t, lower panel) magnitude of the 
crystal field splitting. Results for both "1 down, 2 up" 
(negative Ai) and "1 up, 2 down" (positive Ai) crystal 
field splittings are shown. The chemical potential range 
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FIG. 7: /^-dependence of the orbital occupancy per spin n a (/i) 
for U/t — 12, J/t = 1 at fit = 50 and in the presence of a 
cubic-tetragonal crystal field Ai = ±0.25t (top) and Ai = ±t 
(bottom). The crystal field splits the threefold degenerate d 
level into a doublet and a singlet, with the doublet lying higher 
or lower according to the sign of Ai. The singlet orbital is 
labeled as "orbital 1" and is denoted by open circles (red on- 
line) for Ai < or diamonds (blue on line) for Ai > 0; the 
doublet orbitals are labeled as orbitals 2,3 and are denoted 
by stars (red on-line) for A < or squares (blue on line) for 
A > 0. Magnetic ordering was suppressed by averaging the 
Green function over spin and additional ordering of orbitals 
2 and 3 was suppressed by averaging the Green functions in 
orbitals 2 and 3. The chemical potential range runs from the 
n = 1 Mott phase (fj, ~ 6t) to the n = 3 Mott phase (/i ~ 20t). 
Insulating phases are visible as plateaux in all three densities 
and occur only at integer total density. Orbitally selective 
Mott phases are visible as plateaux in one density with the 
other(s) varying with [i. 



covers dopings from the n = 1 to the n — 3 Mott insulat- 
ing state. The interaction parameters J = t and U = 12t 
are such that the model is insulating at all three of the 
integer fillings n = 1,2,3. 

To discuss the figure it is convenient to begin with 
the "1 down, 2 up", Ai < case (circles and stars, red 
on-line) and to discuss the behavior as the n = 1 Mott 
insulating state found at fi ~ 6t is doped. Although the 
orbital susceptibility of the n = 1 Mott insulating state is 
finite, even the weaker of the two crystal field splittings 
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shown here is larger than the "orbital superexchange" 
and leads to complete orbital polarization. The favored 
orbital is fully occupied (density n — 0.5/spin) and the 
disfavored orbitals are empty. Now consider adding elec- 
trons to the "1 down, 2 up" state. In the weak crystal 
field case (upper panel) we see that (within our resolu- 
tion) the doping-driven Mott transition out of the n = 1 
state is first order: the orbital polarization drops dra- 
matically on doping so that in addition to adding elec- 
trons, doping leads to a transfer of electrons from the 
highly occupied to the less highly occupied orbital. The 
resulting "orbitally polarized Fermi liquid" state evolves 
smoothly upon doping to the obvious two-electron Mott 
state, characterized by the expected 1/2,1/4,1/4 occu- 
pancy per spin. As electrons are added to this two elec- 
tron state, we find a small orbitally selective Mott region 
with band 1 still insulating but bands 2 and 3 metallic. 
At larger chemical potential an insulator-metal transi- 
tion takes place in band 1, leading to an initial decrease 
of the orbital polarization. This state evolves smoothly 
to the three electron, orbitally symmetric state. Thus 
for small crystal field splitting, "orbitally selective Mott 
behavior" only occurs very close to the insulating con- 
centrations. Consider now the larger crystal field split- 
ting (lower panel). In this case the doped state is in 
the orbitally selective Mott phase: at all chemical po- 
tentials the orbital favored by the crystal field splitting 
remains at the Mott occupancy of n = 1 and the carrier 
density varies only in the disfavored orbital, so that one 
has effectively a model of two bands of carriers coupled 
to a spin-1/2 arising from the filled orbital. Similar ef- 
fects were also noticed recently by LiebscbJi in a study 
of Lai-^Sr^TiOa that corresponds to our model in the 
range (0 < n < 1). 

We next turn to the "1 up, 2 down" Ai > case 
(squares and diamonds, blue on-line). At n — 1 we 
see again that the crystal field splitting is large enough 
to fully polarize the Mott insulator (0,1/4,1/4 orbital 
occupancy). In this case, at weak crystal field split- 
ting, the occupancy of the majority orbitals increases 
smoothly with doping (almost all dopants go into the 
initially empty band, leading to a jump at the metal- 
insulator transition). At larger chemical potential, there 
is an abrupt transition to the n — 2 Mott phase with 
0, 1/2, 1/2 orbital occupancy. On further doping to the 
range 2 < n < 3 we observe phenomena analogous to 
those found on doping the n = 1 "1 down, 2 up" state: 
doping leads to a charge transfer between orbitals which 
reduces the degree of orbital disproportionation. At the 
larger crystal field splitting the behavior is different. Be- 
tween n = 1 and n = 2 the minority orbital remains 
empty; the crystal field splitting is large enough to make 
the material effectively a two orbital band insulator. Be- 
tween n — 2 and n — 3 the state is an orbitally selec- 
tive Mott state, with the minority band partially occu- 
pied and coupled to the spin-1 formed by the majority 
states. This physics has also been discussed very recently 
in Ref. [H 
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FIG. 8: (color online) Orbital filling as function of crystal 
field splitting computed for the two electron state with U/t — 
8, pt — 50 and indicated values of J/t. In order to display all 
of the curves on the same figure the crystal field coordinate is 
chosen to be \i — Ai. The red lines with circles correspond to 
the occupancy of orbital 1 and the blue lines with stars to the 
occupancy of orbitals 2, 3. While a metal-insulator transition 
is evident in the curves for J/t = and those for J/t — 
1, Ai > 0, the J/t = 1, Ai < curves exhibit a transition 
to an orbital selective Mott state (band 1 insulating, bands 2 
and 3 metallic). 



VI. CRYSTAL FIELDS IN THE METALLIC 
STATE 

This section considers the effect of crystal field split- 
ting for weaker interactions U < U C 2 where at A Q = 
the system is in the metallic phase. In Fig. [8] we plot 
the variation in orbital occupancies as the crystal field is 
varied at fixed fi corresponding to n — 2, as was done 
in Fig. [5] for a stronger U. The figure shows results ob- 
tained by averaging the Green's functions of orbitals 2 
and 3. For the U = 8t considered here, the Ai = 
metallic phase is characterized by an orbital susceptibil- 
ity Xorb = ~ d ("i-0^ +n3 )/ 2 ) w ith some J dependence but 
a typical magnitude of ~ 0.2 — 0.3 /t. As Ai is increased 
the disfavored orbital 1 empties out and the occupancy 
of the favored orbitals increases. At J = (right hand 
side of the figure) we see that in the "1 down, 2 up" case, 
an apparently first order transition to a (1,0,0) insulating 
state occurs as the magnitude of the crystal field splitting 
increases, whereas in the "1 up, 2 down" case a transition 
occurs to Mott state with two electrons in two orbitals. 
At J = all possible ways of arranging the two electrons 
among the two orbitals are degenerate; the degeneracy 
would be lifted by intersite effects. 

In the more physically relevant J > case (left hand 
side of figure) a crystal field splitting of the "1 down, 2 
up" type leads to an orbitally selective Mott state. On 
the other hand, increasing the magnitude of a crystal 
field splitting of the "1 up, 2 down" type induces a tran- 
sition to a (0, 1/2, 1/2) insulating state. Again, the com- 
putations presented here are for an orbitally disordered 
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FIG. 9: Average energy E of the 2/6 rilled state at U/t = 8, 
J = U/6 as function of crystal field splitting Ai measured rel- 
ative to the energy of the Ai = state. To compensate the 
asymmetry produced by raising/lowering one orbital while 
leaving the other two in place, we have subtracted 2Ai/3£. 
The results with Ai/t < —0.5 arc in an orbital selective 
Mott state with the lower orbital insulating and the other 
two metallic, those in the range —0.5 < A/t < 0.5 are metal- 
lic in all bands and the solutions at > 0.5 are insulating. In 
these phases we observe a linear behavior with slopes 1/3 and 
-2/3. 



Mott state. In particular the orbital selective Mott phase 
would be susceptible to orbital ordering. 

Figure [5] shows the energetics of the orbitally selective 
Mott transition at U = 8t and J « 1.2t with density 
n = 2. The chemical potential has been adjusted to 
keep the occupancy fixed. We see that in the metallic 
phase the energy is hardly affected, while in the orbitally 
selective and insulating phases the energy drops linearly 
with A, with a coefficient given by the occupancy of the 
filled orbital. These energetics are important because in 
several materials (including for example Ca2Ru04) the 
metal insulator transition is of the orbitally selective type 
and is accompanied by a lattice distortion which acts to 
increase the crystal field splitting. 



VII. COMPARISON TO EXPERIMENT 

In this section we place a few relevant materials on 
our calculated phase diagram and discuss implications 
of our results. We begin with SrVC>3, a pseudocubic 
material characterized by 1 electron in the ti g shell, a 
bandwidth corresponding to t w 0.7eV, a U w 5eV w 7t 
and J ~ U/7& SrVC>3 is a good metal, with a modest 
mass enhancement; it is not believed to be close to the 
Mott transition. This behavior is consistent with our 
phase diagram: the value of J/U is similar to that used 
to construct the lower panel of Fig. [2] and U = 7t is quite 
far from the n — 1 Mott lobe. The metallic behavior is 
seen to be a consequence of the non-vanishing value of 
the Hund coupling J. Without J, the material would be 



very close to the Mott transition. In the related material 
LaVO ^ 25 ! 26 ' 27 the change Sr — ► La implies that the d shell 
filing changes from 1 to 2. The material also exhibits a 
moderate orthorhombic distortion away from cubic sym- 
metry, of the "1 down, 2 up" type. Using the U values 
obtained from the singly-occupied system we find that 
the materials would be metallic (albeit in the spin-frozen 
phase discussed in Ref . [H) . 

We next consider LaTiOa, in which the bandwidth is 
such that our parameter t ~ 0.5 — 0.7eV and U ~ 4 — 5eV 
~ 5 — lQt, with J ~ U/6. Here examination of the phase 
diagram reveals that within the single site dynamical 
mean field theory, and in the absence of orbital order- 
ing, the material is not predicted to be a Mott insulator. 
However, it is now known that in the material a substan- 
tial local trigonal distortion occurs i 29 i 30 Our results lend 
support to the idea, advanced in previous papersj&^L that 
the trigonal distortion is essential to the insulating behav- 
ior. The trigonal distortion, by lowering one orbital, will 
effectively convert the problem into a one orbital model. 
Figure 2] shows that the amplitude of the distortion must 
be large, providing a level splitting of the order of one 
quarter of the bandwidth. One difficulty with this in- 
terpretation is that the insulating gap in LaTiC>3 is only 
about 0.2eV, whereas in a single-orbital Mott insulator 
the gap is of order U — 2t ~ 2eV. The small value of the 
gap suggests that effects beyond the scope of the single 
site dynamical mean field theory are important in this 
material. 

We now to turn to the Sr/Ca ruthenates. These are 
materials with two holes in the t^g shell; after a particle- 
hole transformation they map on to the two electron 
case of the model studied here. The "113" materials 
(Si'i_ x Ca 2; )Ru03 crystallize in an orthorhombic struc- 
ture slightly distorted from the cubic perovskite form. 
Both the Sr and Ca materials are metallic, with the Sr 
being ferromagnetic below T c w 150K. The t2 S bands 
have a bandwidth of approximately 2.5eV— correspond- 
ing to a t parameter of about O.QeV in the notations of the 
present paper. The U and J values are not established 
for these compounds but must be substantially less than 
the U = lOt ~ QeV needed to drive a Mott transition. 

The ruthenates also exist in a layered form 
(Sr/Ca)2RuC>4. Here the electron counting is the same as 
in the "113" ruthenates but the tetragonal crystal struc- 
ture means that two of the orbitals have an effectively one 
dimensional dispersion and the remaining one has a two 
dimensional dispersion. Thus, a substantial anisotropy is 
already present in the band structure even in the absence 
of explicit ligand field splitting. At room temperature all 
members of the Sr/Ca series are metallic, but as tem- 
perature is reduced the Ca material undergoes a tran- 
sition to an insulating phase, accompanied by a lattice 
distortion^ while Sr 2 Ru04 remains metallic to lowest 
temperatures. In a recent Letter^ Liebsch and Ishida 
have (in agreement with the prior proposals of Jung et 
al.)2£ argued that the insulating phase should be under- 
stood as a consequence of a "1 up, 2 down" crystal field. 
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Figure \5\ shows that a moderate crystal field splitting of 
the "1-up, 2-down" type can substantially decrease the 
critical U required to drive a metal-insulator transition 
at n = 2. In this interpretation, the effective crystal field 
is small in Sr2Ru04 but increases with Ca-doping, driv- 
ing a metal-insulator transition analogous to that shown 
on the left hand side of Fig. [H The argument in favor of 
a small effective crystal field splitting in Sr2Ru04 is the 
near equality of the occupancies of the d xy , d xz and d yz 
orbitals. A previous weak coupling, Hartree-Fock anal- 
ysis by one of us and Okamoto^ found that at small U 
a non-vanishing J stabilized the systems against crystal 
field distortions. On the basis of this calculation it was 
argued that the near degeneracy of orbital occupancies 
in Sr2Ru04 despite the highly anisotropic crystal struc- 

and the results of 
35l do not survive 



ture was an interaction effect. Figure |5 
Ref. [34| suggest that the results of Ref. 
beyond the weak coupling limit, so that the near equality 
of orbital occupations in Sr2Ru04 must be regarded as 
accidental, with the asymmetry of the bandwidths and of 
the crystal field levels compensating one another to leave 
a small effective splitting. If the near vanishing of effec- 
tive crystal field splitting in Sr2Ru04 is accidental then 
it is very reasonable that the changes induced by Ca sub- 
stitution can move the system away from the accidental 
degeneracy. 

An issue with this interpretation is that at higher tem- 
peratures T > 35QK Ca2Ru04 is metallic, while the stan- 
dard single-site dynamical mean field theory predicts that 
the material should become more insulating as the tem- 
perature is raised. We suggest, following Jung et ai, that 
the temperature-driven first order transition can be un- 
derstood in terms of a temperature-driven lattice distor- 
tion. Indeed the energy diagram (Fig. [9|) in combination 
with a reasonable free energy for lattice distortions, im- 
plies a first order transition. From this figure we see that 
the energy gain per orbital AE from a "1 up, 2 down" 
distortion produced by a crystal field A is approximately 
AE = A — A for A > A with A a fraction of the 
bandwidth parameter t. Let us suppose that a lattice 
distortion of the type observed by Ref. 33 produces a 
crystal field splitting A and that the free energy cost of 
this distortion is AF = A 2 /(2A(T)) w ith A(T) a tem- 
perature dependent free energy scale which increases as T 
is decreased, representing the entropy of thermal phonon 
fluctuations, which favor the undistorted state. The total 
free energy is then 



F = - (A - Ac) 6 (A - A 



1 A 2 



2 A(T) 



(10) 



We see that for A(T) < Ao, F is minimized at A = 
but for A > Ao the free energy is minimized at a value 



A > Ao implying a Mott state, and that as A continues 
to increase the magnitude of A and therefore the size 
of the Mott gap further increases. Additional study of 
this issue using the realistic band structure and a better 
model of the phonon energetics would be desirable. 

Another material to which the considerations of this 
paper should apply is doped Cqq- For this material, 
quantum chemical calculations suggest a E/ ~ 1.5eV, 
J « O.leV^ai and a bandwidth of about 0.6eV^- cor- 
responding in the language of this paper to t ~ 0.15eV 
so U ~ lOt and J ~ {7/15. Experimentally, AiCeo and 
A2C60 are insulators, while A3C60 is metallic and super- 
conducting. The small J limit of the theory is roughly 
consistent with the trend in behavior, with the interac- 
tions being strong enough to place the n = 1 compound 
firmly within the Mott phase while the location of the 
n = 2, 3 materials is uncertain. The calculations pre- 
sented here would suggest that A2C60 and A3C60 should 
be approximately equally likely to be insulating. From 
this point of view, understanding the metallic nature of 
the A3C60 and insulating nature of A2C60 is an impor- 
tant open problem. 



VIII. SUMMARY 

In this paper we have used continuous-time quantum 
Monte Carlo methods to produce a comprehensive pic- 
ture of the metal-insulator phase diagram and response 
to crystal fields of a "three orbital" model which contains 
the essential physics of the fullerides, and the perovskite- 
based titanates, vanadates and ruthenates. We have doc- 
umented the strong effect of the Hunds coupling on the 
location of the Mott transition and on the response to 
crystal fields and have placed a number of experimentally 
interesting materials on the phase diagram. The meth- 
ods presented here provide a basis for detailed, material- 
specific calculations of realistic Hamiltonians. An im- 
portant future direction for research is the investigation 
of the stability of the phases we have found against or- 
bital and magnetic ordering. Work in this direction is in 
progress. 
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